TITLE - Method of Paper Machine Control and Apparatus for the Method 
BACKGROUND OF THE INVENTION 
Field of Invention 

The present invention relates to a method and a system for controlling 
a paper machine, wherein a dryer is controlled by predicting the moisture 
percentage of a web at a dryer part inlet and also predicting the dryer's 
steam pressure according to the predicted moisture percentage. 
Description of Prior Art 

FIG. 1 is a schematic view showing the configuration of a typical paper 
machine. In the figure, raw pulp is discharged from a stock inlet 81 to 
a wire part 82. The wire part 82 is moved in the direction of arrow A 
by means of rotating rolls 821. The raw pulp discharged onto the wire part 
82 is subjected to drainage so as to form a web (that is paper). The web 
thus formed is transferred to a press part 83 for further water drainage. 

The web subjected to water drainage at the press part 83 is transferred 
to a pre-dryer 84. A multitude of steam drums 841 are disposed in the pre- 
dryer 84 and heated by steam introduced thereinto. The web is wound around 
the steam drums as it is moved foward, so that the web is drived until a given 
moisture percentage is reached. 

The dried web is subjected to a sizing process, such as application of 
a sizing agent (coating agent) at a size press 85; is further dried by an 
after-dryer 86; and is then take up as a product indicated by numeral 



87. It should be noted that the after-dryer 86 is configured in the 
same way as the pre-dryer 84. 

Numerals 88 and 89 denote BM systems, both of which detect the basis 
weight, moisture percentage, and other data items of the web as it comes 
out of the pre-dryer 84 and after-dryer 86, respectively. The values of 
data items thus detected are input to a control apparatus not shown in 
the figure. The control apparatus controls the amount of raw pulp 
discharged onto the wire part 82 or the amount of steam introduced into 
the steam drums of the pre-dryer 84 and after-dryer 86, as well as the 
machine speed and other parameters, so that the product in question 
complies with predetermined specifications . Grade change control 
whereby different types of product are produced is also practiced 
commonly . 

In grade change control, any product obtained during the time of 
grade change, wherein a switch is made to another type of product, will 
be treated as broke, i.e., non-standard paper. Therefore, the duration 
of grade change must be minimized in order to increase operation 
efficiency. To solve this problem, an invention of a method of 
predicting a steam pressure setpoint after grade change by simulation is 
described in the specification of Patent 3094798. 

Now, the aforementioned invention is described briefly. 

The invention described in the specification of Patent 3094798 uses 
an iron model wherein the steam drums of the pre-dryer 84 and after- 



dryer 86 are simplified into a planar form. In the model, the state of 
contact among the steam drum, web. and canvas wound continuously round 
the steam drums is classified into five patterns. Then, the heat- 
transfer differential equation of each pattern is derived and converted 
to a difference equation, so that a steam pressure setpoint after grade 
change is predicted by solving the difference equation. 

The heat-transfer differential equations of a pattern wherein the 
steam drum, web and canvas are in contact with each other in this order 
are represented as equations 5 to 7 below. 



1 L »-Po'C D ^ = h s -(T s (t)-T l {t))-h DW -(tM-tM (5) 



K - Pw ■C w d ?M = h DW \T x {f)-TM-Kc iT 2 {t)-TM-Evapo{T 2 ,T w ) -(6) 
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The meanings of the parameters included in equations 5 to 7 are as 
follows . 

L 0 : Drum thickness (m) 
L v : Web thickness (m) 



L c : Canvas thickness (m) 

T,: Steam temperature. within drum (°C) 

T,: Drum's surface temperature (°C) 

T 2 : Web (paper) temperature (°C) 

T 3 : Canvas temperature (°C) 

T a : Dry-bulb temperature of air within hood (°C) 

C D : Drum's specific heat (kJ/(kg°C)) 

C„: Web's (paper's) specific heat (kJ/(kg-°C)) 

P 

jS C c : Canvas' specif ic heat (kJ/(kg-°C)) 

m 

J; p D : Drum's density (kg/m 3 ) 
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p w : Web's (paper's) density (kg/m ) 
p c : Canvas' density (kg/m 3 ) 



Hi h s : Coefficient of heat transfer between steam within drum and drum 

DO ■ 

Pi 2 

y surface (kJ/ (m -sec- C ) ) 

h w : Coefficient of heat transfer between drum surface and web 
(kJ/ (m 2 -sec-°C) ) 

- h^: Coefficient of heat transfer between web surface and canvas 

(kJ/(m 2 -sec-°C) ) 

h a : Coefficient of heat transfer between canvas and air within hood 
(kJ/(ra J -sec-°C) ) 

FIG. 2 is a table that summarizes the above-listed parameters. 

The term Evapo(T 2 , T H ) in equation 6 is a function representing the 
amount of heat of evaporation removed from the web as the result of 
moisture evaporation, and is- given by equation 8 below. 



Evapo(T 2 , TJ = V(MP ABS )-K-(P(T 2 )-P(T W ))-SB(T 2 ) (fc//(m 2 -sec)) (8) 



where 

P(T) = Saturation vapor pressure (kPa) at temperature T (°C) 
SB(T) = Heat of evaporation (kJ/H 2 0kg) at temperature T (°C) 
T M = Wet -bulb temperature of air within hood (°C) 

VtMP^) = Function representing moisture evaporation intensity at 
absolute moisture percentage MP^, where 0.0 <; VfMP^) <; 1.0 
( dimensionless ) 

K = Drying rate coefficient ( H 2 0kg/ (rn 2 -sec-kPa) ) 

Although heat-transfer differential equations for patterns of contact 
other than those mentioned above are also given by the invention 
described in the specification of Patent 3094798, these equations are 
omitted here to avoid complication. 

In differential equations 5 to 7 discussed earlier, a length of time is 
segmented into time intervals At, which is determined by the machine 
speed, circumference of a steam drum, and other data items, so that a 
difference equation is derived and the numeric solution thereof is 
obtained. Since the web moves from the upstream side to the downstream 
side of the paper machine as time elapses, it is possible to calculate 
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the web temperature at the steam drum by numerically solving the 
difference equation. 

From equation 8, EvapoMP(T 2 , T H ) (H 2 Okg/ (m 2 -sec) ) , which is the amount 
of moisture evaporated from the web per unit area and unit time, can be 
represented by equation 9 below. 



By using this equation, it is possible to calculate the absolute 



of the incremental time interval At as shown in equation 10 below. 



Eva P oMP(T 2 , tj = v(mP abs )-K-(p(t 2 )-P{T w )) [h 2 Okg l(m 2 sec] 



(9) 



moisture percentage MP^sCj) (j = 1 



N) of the web after the lapse 



MP ABS (j * U = MP ABS (j). 



10 3 • EvapoMP(T 2 ,T w ) • At 
BD 



(10) 



where 



BD = Bone-dry basis weight (g/m 2 ) 



At = 



Incremental time interval (sec) 



MP* BS (j) (3 = 1 



N) = Absolute moisture percentage (%) at mesh 



division j 
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From this absolute moisture percentage, it is possible to calculate 



MP(j) (j = l N) = Relative moisture percentage (%) at mesh 

division j 

FIG. 3 is a flowchart representing the algorithm of a steady- state 
simulation using equations 5 to 11. In the first step, the current 
operation status data, such as the current machine speed (m/min), basis 
weight setpoint (g/m 2 ), and moisture percentage setpoint (%), are 
acquired. In the second step, the incremental time interval At for 
differential calculations is determined from the machine speed, drum's 
circumference, and other data items. In the third step, the steam 

temperature T s (j) (j = 1 N) within the drum is calculated from 

the current dryer steam pressure setpoint by using a saturation vapor 
pressure curve. Note that N is the number of mesh divisions. 

In a further step, equations 5 to 11 and the difference equations 
derived therefrom are used to calculate the drum temperature T t (j) (j = 
1, . . , N) , web temperature T 2 (j) (j = 1, . . , N), canvas temperature 



the (relative) moisture percentage MP(j) (j = 1 



N) (%) as shown 



in equation 11 below. 




(11) 



where 
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T 3 (j) (j = 1, . . , N), and web's final moisture percentage MP(j) (j = 
1, . - , N). In yet a further step, a judgment is made on convergence 
between the web's relative moisture percentage MP(N) and actual measured 
value MP MEASURE provided by a moisture sensor at a final cylinder. 
Convergence has been reached if the absolute value of the difference 
between MP(N) and MP^^g is smaller than the given value EP. 

If convergence has not yet been reached, the drying rate coefficient 

jMi K is corrected by AK to calculate the drum temperature, web temperature, 

D 

p canvas temperature, and web's relative moisture percentage once again. 

yl 

WjJ When convergence has been reached, the drying rate coefficient K, drum 

m 

y ; temperature Tj^j), web temperature T 2 (j), canvas temperature T 3 (j), and 
O web's moisture percentage MP(j) are fixed to their values at that moment, 

y, 

W and the steady-state simulation ends. 

Si 
| 

fU For a dryer part consisting of pre-dryer and after-dryer parts, it is 

also acceptable to calculate the moisture percentage at an after-dryer 
outlet as the final moisture percentage. Alternatively, moisture 
percentages at the pre-dryer and after-dryer outlets may be defined as 
the final moisture percentages. In the latter case, a convergence 
calculation should be made for each of the dryer parts . 

In the steady-state simulation heretofore discussed, the drying rate 
coefficient K is adjusted so that the absolute moisture percentage at 
the final cylinder is approximated to the actual measured value. Next, 
a simulation of steam pressure prediction is carried out, in order to 



predict the optimum steam pressure setpoint in an operation status after 
grade change. The simulation of steam pressure prediction is explained 
by referring to FIG. 4. 

In the first step in FIG. 4, operation status data after grade change, 
i.e., the machine speed (m/min), basis weight setpoint (g/m 2 ), and 
moisture percentage setpoint (%), are acquired. In the second step, the 
incremental time interval At for differential calculations is determined 
from the machine speed, drum's circumference, and other data items. In 

the third step, the steam temperature T s (j) (j = 1 N) within the 

drum is calculated from the current dryer steam pressure setpoint P 
(kPa) by using a saturation vapor pressure curve. Note that N is the 
number of mesh divisions . 

In a further step, the value of the drying rate coefficient K 
determined in the steady-state simulation, as well as the value before 
grade change used in the steady-state simulation, for example, as the 
pre-dryer part inlet moisture percentage, is used to find the numerical 
solutions of equations 5 to 11 and their difference equations, thereby 
calculating the drum temperature T x (j) (j = 1, . . , N), web temperature 
T 2 (j) (j = 1, - • , N) , canvas temperature T 3 ( j ) (j = 1, . . , N), and 
web's moisture percentage MP(j) (j = 1, . . , N) as the initial values 
for the difference equations . 

In yet a further step, the value of the web's moisture percentage 
MP(N) at the final cylinder and the moisture percentage setpoint after 



grade change are compared, in order to judge convergence in the same way 
as in the case of the steady- state simulation. If convergence has not 
yet been reached, the dryer steam pressure setpoint is corrected by the 
given value At, and the drum temperature, web temperature, canvas 
temperature, and web's relative moisture percentage are calculated once 
again. When convergence has been reached, the values of these data 
items at that moment are fixed and the simulation of steam pressure 
prediction ends. 

In such a paper machine as discussed above, controlling the process 
of drying a product is an important factor in order to produce products 
of consistent quality. Drying at the after-dryer 86 is particularly 
important since the drying process directly affects product quality. 
For this reason, it is necessary to precisely know the moisture 
percentage of a product at the dryer inlet. 

Traditionally, the moisture percentage of a product at the inlet of 
the after-dryer 86 has been calculated by using a measured value 
provided by the BM system 88 installed before the size press 85 and then 
applying, for example, equation 12 shown below. It should be noted that 
the absolute moisture percentage in the equation means the ratio of 
moisture weight to the bone-dry weight of a web which is a product. 



BD. 



x absMP p 



+ CW- 



100-5 



PRE 



PREEND 



s_ 



absMP A 



(12) 



BD f 



'aft 



where 

absMPj^m = Absolute moisture percentage (0.0 to 1.0) at after-dryer 86 
inlet 

absMP PREEND = Absolute moisture percentage (0.0 to 1.0) at pre-dryer 84 
outlet (calculated by simulation) 

BD PRE = Bone-dry basis weight (g/m 2 ) at pre-dryer 84 outlet (measured with 
BM system) 

b £>aft = Bone-dry basis weight (g/m 2 ) at after-dryer 86 outlet (measured 
with BM system) 

Ctf = Size's bone-dry coated weight (g/m 2 ) 

S = Moving average of size's (coating agent's) concentration (%) 

The pre-dryer 84 outlet absolute moisture percentage absMP PREEND is 
evaluated as a solution given by simulating the steady state formed in 
the pre-dryer 84. However, a size with a concentration of 5 to 10% is 
coated at the size press 85 and therefore, the moisture percentage must 
be corrected by the amount of moisture produced by such coating. 

More specifically, the first term BD PRE x absMP PREEND of the numerator 
on the right-hand side of equation 12 denotes a moisture weight (g/m 2 ) 
per unit area at the outlet of the pre-dryer 84, whereas the second term 
CW-(100 - S)/S denotes a moisture weight (g/m 2 ) contained in the coated 
size per unit area. Since the sum of these two terms is the amount of 
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moisture contained per unit area of a product at the inlet of the after- 
dryer 86, it is clear that the absolute moisture percentage is evaluated 
by dividing this amount by the bone-dry basis weight BD^ measured with 
the BM system 89. 

It should be noted that as the size's bone-dry coated weight CW, 
equation 12 uses the value calculated by equation 13 below, which is the 
difference between the bone-dry basis weights measured with the BM 
systems 88 and 89. 

CW = BDjpp - BD PRE 

(13) 

SUMMARY OF THE INVENTION 

The following problems have been inherent, however, with the method 
of calculating the moisture percentage of a web at a dryer inlet in such 
a paper machine as discussed above and with the simulation of steam 
pressure prediction after grade change. 

In the simulation of steam pressure prediction shown in FIG. 4, the 
prior art method uses, for example, the same moisture percentage as that 
before grade change, i.e., a value input from a moisture sensor or such 
a fixed value as 50%, for the initial moisture percentage MP(1). If 
operation status data, such as the basis weight or machine speed, 
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changes before or after grade change, then there will also be a change 
in the wire retention which is the ratio at which raw material remains 
on a wire part, in the concentration of circulating white water, or in 
the capability of water drainage at a wire part. Accordingly, it is 
known that the moisture percentage (normally 40 to 60%) at the dryer 
inlet also changes consequently. 

Furthermore, it is empirically known that the moisture percentage at 
the dryer inlet increases if the basis weight increases while the 
machine speed is kept constant. For example, if the basis weight 
changes by 10 g/m 2 , the moisture percentage changes by 1 to 2% . If a 1% 
change takes place in the moisture percentage at the dryer inlet , there 
will be an approximately 10 kPa change in the predicted value in the 
simulation of steam pressure prediction shown in FIG. 9. For this 
reason, the prior art method has the problem that if the moisture 
percentage at the dryer inlet is set to the same value as that before 
grade change, the predicted value of steam pressure contains a non- 
negligible error, against a desirable steam pressure setpoint after 
grade change. 

As is clear from equation 12, the bone-dry basis weight BD PRE at the 
pre-dryer 84 outlet is used to calculate the moisture percentage. 
Consequently, the prior art method has another problem that in the case 
of a paper machine not provided with the BM system 88 that should 
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otherwise be located before the size press 85, it is not possible to 
calculate the moisture percentage using the equation. 

Yet another problem with the prior art method is that even if the BM 
system 88 is installed at all, the error in moisture percentage 
calculation due to the measurement accuracy of the BM system tends to 
become large if the moisture percentage is calculated using equations 12 
and 13. In other words, the method of calculation using equations 12 
and 13 is problematic since it involves significant calculation errors. 
This problem is explained here by citing specific examples . 

Now, we assume that individual measured values and the moisture 
percentage calculated therefrom are as follows. 

BD PRE = 100.0 (g/m 2 ) 

BD AFT = 102.0 (g/m 2 ) 

CW = 2.0 (g/m 2 ) 

S = 8% 

absMP PREEND =0.02 

By substituting these values into equation 12, we obtain 
BD PRE x absMP PRESND = 100 x 0.02 = 2.0 
CW-(100 - S)/S = 2 x 11.5 = 23.0 
absMP^^ = (23.0 +• 2.0)/102.0 = 0.245 



On the other hand, the accuracy ranges of individual measuring 
instruments are approximately as follows . 



Accuracy range of basis weight sensor: ±0.15 (g/m 2 ) 



Accuracy range of moisture sensor: ±0.1 (%) 



From these values, the accuracy levels of bone-dry basis weight and 
bone-dry coated weight can be calculated as shown below. 



b 

[J Accuracy of bone-dry basis weight =V0.1x 0.1 + 0.15 x 0.15 = 0.18 
1 " ™ 

■ jj Accuracy of bone-dry coated weight ACW =^0.18 X 0.18 + 0.18 x 0.18 = 0.25 

u 

g«j From these calculations, errors in the size's coated weight per unit 

m area and in the moisture percentage at the after-dryer 86 inlet are 

C determined as shown below. 

jl) 



Accuracy of size's coated weight ACW •— — — = 0.25x11.5 = 2.88 

S 

Accuracy of moisture percentage AabsMP AFTIN at after- dryer 86 inlet 

isr^— — - °-° 28 

BD AFT 102.0 

This means that an error as large as AabsMP AFTIN /absMP AFT i N = 0.028/0.245 
= 11.4% will occur depending on the accuracy of measuring instruments. 
The error in the bone-dry coated weight is comparatively large since the 
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weight is calculated from measured values provided by four measuring 
instruments. This error increases further by a factor defined by the 
reciprocal of the size's concentration (approximately 10 times larger), 
if the amount of moisture contained in the size is included in the 
calculation. Consequently, such a large error as mentioned earlier 
occurs . This error results in yet another problem that precise control 
is not feasible. 

As discussed heretofore, it is evident that precisely predicting the 
pi ; moisture percentage at a dryer inlet is of great significance in paper 
igl machine control. 

if! 

M. It is therefore an object of the present invention to provide a 

P method of paper machine control whereby the moisture percentage of a web 

U 

h i 

lz at a dryer inlet is estimated, excellent control can be performed, and 

W 
O 

pi the time required for grade change can be reduced, as well as apparatus 
for the method. 



BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic view showing the configuration of a typical 
paper machine. 

FIG. 2 is a table that summarizes parameters used in heat- transfer 
equations . 
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FIG. 3 is a flowchart representing steady-state simulation in a prior 
art method. 

FIG. 4 is a flowchart representing the simulation of steam pressure 
prediction in the prior art method. 

FIG. 5 is a flowchart representing one embodiment in accordance with 
the present invention. 

FIG. 6 is a flowchart representing another embodiment in accordance 
with the present invention. 

FIG. 7 is a block diagram showing the configuration of one embodiment 
in accordance with the present invention. 

FIG. 8 is a block diagram showing the configuration of yet another 
embodiment in accordance with the present invention. 

FIG. 9 is a block diagram showing the configuration of yet another 
embodiment in accordance with the present invention. 

FIG. 10 is a diagrammatic view showing the configuration of yet 
another embodiment in accordance with the present invention . 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Preferred embodiments of the present invention will now be described 
in detail with reference to the accompanying drawings. 



A method of evaluating a steam pressure setpoint after grade change 
at a'pre-dryer part is first explained by referring to FIGS. 3 and 5. 
In the first step, steady-state simulation is performed in the same way 
as for the prior art method shown in FIG. 3, in order to determine the 
drying rate coefficient K. In the second step, steam pressure after 
grade change is predicted. As discussed earlier, a value before grade 
change (e.g., 50%) is used as the moisture percentage MP(1) at a dryer 
part inlet in the prior art method of simulating the prediction of steam 
pressure after grade change. This method of prediction results in the 
problem that an error occurs in the predicted value of steam pressure. 
For this reason, equation 18 shown below is used to calculate the 
moisture percentage at the dryer part inlet (pre-dryer part inlet), 
which is the web ' s initial relative moisture percentage MP ( 1 ) , when . 
making numeric calculations for solving difference equations based on 
heat-transfer differential equations in the simulation of steam pressure 
prediction represented by the flowchart of FIG. 5. 

MPNowInit + 4 • BL>2 ~ BD > + A 2 • (18) 

where 

MPNowInit = Initial moisture percentage (e.g., fixed to 50%) at dryer 
part inlet 



MPNextlnit = Initial moisture percentage at dryer part inlet for 
simulation of steam pressure prediction 

BD 2 = Bone-dry coated weight (g/m 2 ) before grade change 

BD 2 = Bone-dry coated weight setpoint (g/m 2 ) after grade change 

V z = Machine speed (m/min) before grade change 

V 2 = Machine speed setpoint (m/min) after grade change 

Aj = Ratio of change in dryer inlet moisture percentage to change in 
basis weight 

A 2 = Ratio of change in dryer inlet moisture percentage to change in 
machine speed 

MPNowInit is also the initial value of the dryer part inlet moisture 
percentage used in the steady-state simulation shown in FIG. 3. The 
measured value of MPNowInit may be used if a moisture sensor is 
installed. It is also acceptable to use a fixed value or any other 
value appropriate for the condition of operation 'in the absence of the 
moisture sensor. For example, a value (e.g., '50%) input from the 
moisture sensor may be used. Note that A x and A 2 are tuning parameters 
and are adjusted according to the operation status. If we assume that 
A : = A 2 = 0, the result of the simulation equals that of the prior art 
method of simulation shown in FIG. 4. Also note that MPNowInit can 
likewise be used as a tuning parameter. 



Now, an example of calculation based on equation 18 is shown. If we 
define A x and A 2 as A l = 8.7 and A 2 = 1.0 and assume that BD L = 100 (g/m 2 ) , 
BD 2 = 125 (g/m 2 ), Vj = 750 (m/min), and V 2 = 600 (m/min), then equation 18 
is calculated as 

MPNextlnit = 50 + 8.7 x 25/100 + 1.0 x (-150)/750 = 52.0% 

This calculation example shows that there is a 2% difference from the 
case wherein the initial value of 50% before grade change is used as an 
initial value after grade change for steady-state simulation. 

FIG. 5 is a flowchart representing the simulation of steam pressure 
prediction after grade change performed using equation 18. 

In FIG. 5, operation status data, such as the machine speed, is read 
first. Then, an incremental time interval is determined from the 
machine speed, circumference of the drum, and other data items. These 
steps are the same as those of the prior art method of FIG. 4. Next, 
MPNextlnit is calculated using equation 18 and the resulting value is 
substituted into the initial dryer inlet moisture percentage MP(1). In 
a further step, the drum temperature, web temperature, canvas 
temperature, and web's moisture percentage are calculated, in order to 
examine whether or not the web's final moisture percentage MP(N) has 
reached convergence. If convergence has not yet been reached, the steam 
pressure setpoint is corrected to once again calculate the drum 
temperature, and so on. This step is also the same as that of the prior 
art method of FIG. 4. The steam pressure setpoint thus corrected is 



used as the dryer steam pressure setpoint after grade change in order to 
control the paper machine. 

FIG. 7 is a block diagram showing the configuration of a paper 
machine control system whereby the method of paper machine control, 
including the method of steam pressure prediction shown in the flowchart 
of FIG. 5, is implemented. In FIG. 7, numeral 11 denotes an initial 
settings block for reading the current operation status or determining 
the incremental time interval At. Numeral 12 denotes a relative 
moisture percentage calculation block for calculating the initial value 
of a moisture percentage according to equation 18 discussed earlier. 
Numeral 13 denotes a drying rate coefficient calculation block for 
evaluating the drying rate coefficient by simulation according to the 
flowchart shown in FIG. 3. 

Numeral 14 is a steam pressure prediction block, to which the outputs 
of the initial settings block 11, relative moisture percentage 
calculation block 12, and drying rate coefficient calculation block 13 
are applied, in order to make a loop calculation in the algorithm 
flowchart of FIG. 5 and predict steam pressure after grade change. 
Numeral 15 is a controller for controlling the paper machine by using 
the steam pressure predicted by the steam pressure prediction block 14 
as a dryer steam pressure setpoint after grade change. Numeral 16 is a 
dryer to be controlled. Thus, it is possible to realize a paper machine 
control system that enables the duration of grade change to be reduced. 
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Now, a method for evaluating a pressure setpoint at an after- dryer 
after grade change is explained. Firstly, steady-state simulation is 
explained by referring to FIG. 3. In the simulation, the absolute 
moisture percentage absMP^™ at an after-dryer inlet in a steady state 
is calculated using equation 12 discussed earlier. At this point, note 
that absMP PREEND is determined by the steady-state simulation of the pre- 
dryer. In addition, the bone-dry coated weight CW of a size is 
determined from equation 19 shown below. 

1000xFx{S/100)xW ) 
Vxd 

where 

CW = Size's bone-dry coated weight (g/m 2 ) 

F = Moving average of size's flow rate (L/min) 

S = Moving average of size's concentration (%) 

W = Size's specific gravity (kg/L) 

V = Machine speed (m/min) 

c? = Paper width (m) 



The numerator of equation 19 is a product of the size's flow rate and 
concentration, thus representing the bone-dry weight of the size 
consumed in one minute. The numerator therefore has a unit of g/minute. 
The concentration S, which has a unit of %, is divided by 100 so that it 
is converted to a ratio. Likewise, the specific gravity W, which has a 
unit of kg/L, is multiplied by 1000 so that the unit is converted to 
grams . 

The denominator of equation 19 is a product of the machine speed and 
paper width, thus representing the area of paper onto which the size is 
transferred in one minute. The denominator therefore has a unit of 
m 2 /minute. Accordingly, by using this equation it is possible to 
determine the weight of the size transferred onto unit area of paper, 
i.e., the size's bone-dry coated weight CW. 

The size's flow rate F and concentration S are measured with a 
flowmeter and concentration meter, respectively. Thus, the moving 
averages of these parameters are taken over a sufficiently long period 
of time such as five minutes since the parameters are not for use in 
quick-response, dynamic control. For this reason, it is possible to 
minimize the effect of short -period variations or errors in the measured 
values of the parameters even if there is any such variation or error. 

From the absolute after-dryer inlet moisture percentage absMP AFTIN , the 
web's initial moisture percentage MP(1) in the calculations for finding 
the numeric solutions of the difference equations in FIG. 3 is 
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calculated according to equation 11. This calculation determines the 
drying rate coefficient K at the after-dryer part, as shown in FIG. 3. 
If a BM system is also installed before a size press, a convergence 
calculation is made separately for the pre-dryer and after-dryer. If 
not, a convergence calculation is made for the after-dryer only. 

FIG. 8 is a block diagram showing the configuration of a paper 
machine control system that uses equation 19. In the figure, numeral 21 
denotes a web production block for producing a web not yet subjected to 
size coating. Numeral 22 denotes a size coating block for coating the 
size onto the web produced at the web production block 21. Numeral 23 
denotes a dryer for drying the web onto which the size has been coated. 
Numeral 24 denotes a moisture percentage calculation block for 
calculating the moisture percentage of the web coated with the size from 
equations 12 and 19 discussed earlier. Numeral 25 denotes a controller, 
to which the moisture percentage calculated and output by the moisture 
percentage calculation block 24 is input in order to control the dryer 
23 according to the moisture percentage. 

Now, the simulation of after-dryer steam pressure prediction during 
grade change shown in FIG. 6 is discussed. Under normal conditions, the 
machine speed changes during grade change. Since the amount of size 
transferred at the size press is proportional to the machine speed, the 
size's flow rate also changes in proportion to the machine speed if the 
machine speed changes. Consequently, equation 20 shown below holds true 



where 

F = Moving average of size's flow rate (L/min) at machine speed V 

F* = Moving average of size's flow rate (L/min) at machine speed V* 

V, V = Machine speed (m/min) 

If the machine speed, changes from V to V* and the size's 
concentration from S to S* through the grade change, then the size's 
bone-dry coated weights CW and CW* before and after the grade change can 
be calculated according to equation 21 below. 



1000 xF* x[S' /lOOjxiy 
V'xd 

1000xFx(S/100)xH/ S' 
Vxd S 

S' 



(21) 



where 

CW and CW = Bone-dry coated weights (g/m 2 ) before and after grade change, 
respectively 
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F and F* = Moving averages of size's flow rates (L/min) before and after 
grade change, respectively 

S and S* = Moving averages of size's concentrations (%) before and after 
grade change, respectively 

W = Size's specific gravity (kg/L) 

d = Paper width (m) 

According to equation 21, it is possible to predict the bone-dry- 
coated weight after grade change using equation 22 below if 
concentration setpoints of the size are given for each grade. 



CW-^- (22) 



where 

CW = Bone-dry coated weight (g/m 2 ) before grade change based on equation 
19 

CPf = Predicted bone-dry coated weight (g/m 2 ) after grade change 

S T and S' T = Size's concentration setpoints (%) before and after grade 
change , respectively 

This means that it is possible to know the bone-dry coated weight 
after grade change before a grade change takes place. As explained 
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earlier, the absolute pre-dryer outlet moisture percentage absMP PEtEEND * 
after grade change can be evaluated by simulation. In addition, the 
pre-dryer outlet bone-dry basis weight is given by subtracting the 
size's bone-dry coated weight from the after-dryer outlet bone-dry basis 
weight. Consequently, it is possible to predict the absolute after- 
dryer inlet moisture percentage absMP AFTIN * after grade change before a 
grade change takes place, according to equation 23 below. 



. 100-5; 
(BD AfT -CW )xabsMP PREEND +CW 



••(23) 



where 

absMP^^' = Absolute moisture percentage (0.0 to 1.0) after grade change 
at after-dryer 86 inlet 

absMP PREBtlD ' = Absolute moisture percentage (0.0 to 1.0) after grade change 
at pre-dryer 84 outlet 

BD AFT = Bone-dry basis weight (g/m 2 ) at after -dryer 86 outlet (measured 
with BM system) 

CfV* = Predicted bone-dry coated weight (g/m 2 ) after grade change 

S' T = Size's (coating agent's) concentration setpoint (%) after grade 
change 
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By referring to FIG. 6, a method for calculating the predicted, value 
of after-dryer steam pressure after grade change using the absolute 
after-dryer moisture percentage absMP^^* after grade change is 
explained. In FIG. 6, the steady-state steam pressure setpoint of each 
section of the dryer before grade change, and operation status data 
after grade change, such as the machine speed, are read first. Then, an 
incremental time interval is determined from the machine speed, 
circumference of the drum, and other data items. These steps are the 
same as those of the prior art method of FIG. 4. In a further step, 
absMP AFTIN * is calculated using equation 23, and the resulting value is 
substituted into the initial dryer inlet moisture percentage MP(1). In 
yet a further step, the drum temperature, web temperature, canvas 
temperature, and web's relative moisture percentage are calculated, in 
order to examine whether or not the web's final moisture percentage 
MP(N) has reached convergence. If convergence has not yet been reached, 
the steam pressure setpoint is corrected to once again calculate the 
drum temperature, and so on. This step is also the same as that of the 
prior art method of FIG. 4. The steam pressure setpoint thus corrected 
is used as the dryer steam pressure setpoint after grade change in order 
to control the paper machine. 

A method of switching from the steam pressure setpoint before grade 
change to the above-mentioned steam pressure setpoint after grade change 
may be in compliance with the method described in the specification of 



Patent 3094798 filed earlier. Other alternative methods may also be 
permissible. 

FIG. 9 is a block diagram showing the configuration of a system for 
paper machine control at the time of grade change. In the figure, 
numeral 31 denotes a web production block for producing a web not yet 
subjected to size coating. Numeral 32 denotes a size coating block for 
coating the size onto the web produced at the web production block 31. 
Numeral 33 denotes a dryer for drying the web onto which the size has 
been coated. Numeral 34 denotes a moisture percentage prediction block 
for predicting the moisture percentage of the web after grade change 
from equations 12 and 22 discussed earlier. Numeral 35 denotes a 
controller, to which the output of the moisture percentage prediction 
block 34 is applied in order to control the dryer 33. The controller 35 
controls the dryer 33 according to the output of the moisture percentage 
prediction block 34 after the grade change takes place. 

FIG. 10 is a diagrammatic view showing the configuration of a system 
for controlling size coating. In the figure, numeral 4 denotes a 
concentration controller, wherein a size with a constant concentration 
stored in a storage tank 5 is mixed with dilution water to produce a 
size with a desired concentration. The flow rate of a size with a 
constant concentration of C (= 10%) is detected by a flowmeter 46 and 
input to a ratio setting unit 41. Also, the ratio of dilution water is 
input to the ratio setting unit 41 manually for each grade. The ratio 
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setting unit 41 controls a valve 42 so that a preset ratio of dilution 
water is reached. The^size with the constant concentration is mixed 
with dilution water at a rotary screen 61 and stored in a supply tank 62. 

The size stored in the supply tank 62 is injected into a coater 71 
within a size press 7, transferred to a roll 72, and then transferred 
further to a web (paper) 73 which is a product. The level of the supply 
tank 62 is measured with a level meter 63 and the measured value is 
input to a valve controller 44. The valve controller 44 controls a 
valve 45 so that the level of the supply tank 62 is kept constant. 

Since the rate of transfer at the coater 71 is constant, the flow 
rate of the size changes in proportion to the machine speed. 

If we assume the flow rate of the size supplied from the storage tank 5 
is A (L/min) and the ratio set in the ratio setting unit 41 is r, then 
the flow rate F (L/min) of the size supplied to the coater 71 is 

F = (1 + r)-A 

In addition, the following relationship exists between the size's 
concentration S and the dilution water's ratio r. 



{l + r)xA 1 + r 
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Consequently, it is possible to calculate the bone-dry coated weight 
CW from equation 19 and evaluate a bone-dry coated weight after grade 
change from equation 22. According to the results of the calculation 
and evaluation and the result of calculating a pre- dryer outlet moisture 
percentage based on the simulation of the condition of drying by the 
pre-dryer, it is also possible to calculate the after-dryer inlet 
moisture percentage before grade change from equation 12 discussed 
earlier. Furthermore, it is possible to evaluate a moisture percentage 
after grade change by substituting the CW* T of equation 22 for the CW of 
equation 12 . 

As is evident from the description heretofore given, the following 
advantages can be expected according to the present invention . 

In one aspect of the paper machine control method according to the 
present invention, wherein a dryer steam pressure after grade change is 
predicted by solving difference equations obtained by differentiating 
heat-transfer equations that hold true among a steam drum, web and 
canvas and the predicted value is used as a dryer steam pressure 
setpoint after grade change, the initial value of a relative moisture 
percentage at a dryer part (pre-dryer part) inlet is calculated 
according to a given equation when solving the difference equations . 

Accordingly, it is possible to obtain a value closer to an actual 
steam pressure setpoint as the predicted value of a dryer steam pressure 
after grade change. Consequently, it is possible to reduce the duration 



of grade change by adopting the predicted value as the steam pressure 
setpoint after grade change,- reduce the amount of broke; and improve 
productivity. 

Another advantage is that since such items of data concerning the 
drying condition within the dryer as the web temperature and moisture 
percentage can be predicted with higher precision, it is possible to 
provide an operator with more useful information for operations. 

In another aspect of the present invention, the parameters A 1( A 2 and 
MPNowInit are tuned according to the operation status. Consequently, 
the versatility of the paper machine control method increases , since the 
method can deal with a variety of paper machines and operation statuses. 
The versatility can further be increased by tuning the parameters. 

In yet another aspect of the present invention, the bone-dry coated 
weight of a size is calculated according to a given equation; the 
moisture percentage of a web at an after- dryer part is predicted from 
the bone-dry coated weight; and the dryer is controlled using the 
predicted moisture percentage. 

Consequently, it is possible to precisely calculate the coated weight 
even if no BM system is installed before the size press. This means 
that the dryer can be controlled easily by measuring, only the moisture 
percentage at the after- dryer part outlet and making a convergence 
calculation. It is also possible to control the dryer with higher 
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precision since a precise coated weight can be evaluated without being 
affected by instrument errors, thereby improving product quality. 

Yet another advantage is that the control method can be used for 
operation monitoring or steady- state control if there are no BM systems 
installed. 

Yet another advantage is that an apparatus for the control method can 
be built more easily and economically if the number of BM systems can be 
reduced . 

Yet another advantage is that it is possible to precisely estimate 
the moisture percentage after grade change, thus reducing the duration 
of grade change and the amount of broke and improving productivity. 

In yet another aspect of the present invention, the moving averages 
of measured values are used as the flow rate and concentration of a size. 
Consequently, it is possible to prevent the effect of short-period 
variations or errors in flowmeters and concentration meters , whereby the 
moisture percentage can be estimated with higher precision. Furthermore, 
it is possible to use inexpensive flowmeters and concentration meters. 
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